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Abstract 

We calculate the structure properties of protoneutron star such as equation of state, maximum 
mass, radius and temperature profile using the lowest order constrained variational method. We 
show that the mass and radius of protoneutron star decrease by decreasing both entropy and 
temperature. For the protoneutron star, it is shown that the temperature is nearly constant in the 
core and drops rapidly near the crust. 
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I. INTRODUCTION 



Neutron star which is highly compact stellar objects is a result of supernova explosion. 
The structure properties of this object, especially its maximum mass, is of a great interest 
for astrophysicists. Since the compactness parameter for a neutron star is about 0.2 — 0.4, its 
structure should be studied using the general theory of relativity. In fact, the computation 
of the structure properties of a neutron star can be derived using the Tolman-Oppenheimer- 
Volkoff (TOV) equation 

Just after the supernova collapse, a newly-born neutron star called protoneutron star is 
formed. At these stages, the neutron stars are rich in leptons. This is due to the fact that the 
neutrinos are trapped in the protoneutron star matter. The temperature of protoneutron 
.a. . g.ate. than lOMeV ^ fl fl B fl B 1 The.efo.. the ten,pe.atu. of ttae 
stages cannot be neglected with respect to the Fermi temperature throughout the calculation 
of its structure. Depending on the total number of nucleons, a protoneutron star evolves 



either to black hole or to stable neutron star 



0,y, 3- 



Therefore, computation of the 



maximum mass of protoneutron star is of crucial importance. 



n recent years, we have investigated the properties of neutron star matter 



3, 



11 



12, 



14| . In the present work, we intend to compute the structure properties of protoneutron 
star at different stages employing the lowest order constrained variational approach using 
the modern AVis potential |15| . 



II. FORMALISM 



In this section, we give the formalism of calculation for the protoneutron star matter 
equation of state. 

We write the total energy per baryon (E) as the sum of contributions from leptons and 
nucleons, 

E = Eiep + Enud ■ (1) 

The contribution from the energy of leptons is 

Elep = „ / n{x)Vl + x^x^dx, (2) 
n'^pn Jo 

where p is the total number density of nucleons and n{x) is the Fermi-Dirac distribution 
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function, where x defined, 

x = (3) 

mc 

We calculate the contribution from the energy of nucleons using the lowest order constrained 
variational method. We consider up to the two-body term in the cluster expansion for the 
energy functional, 

Enud = Ei + E2- (4) 

The one-body energy Ei is 

The two-body energy E2 is 

E2 = {2Ar'Y.< ^3 I {-7^ [/(12), [Vi2,/(12)]] + /(12)V(12)/(12)} | ij >„ (6) 
where f(12) and V(12) are the two-body correlation function and nucleon-nucleon potential. 



respectively. V(12) has the following general form |15| . 

18 

y(12) = ^V>i2)0?2- (7) 
p=l 

Finally, we minimize the two-body energy, i?2; with respect to the variation in the two-body 
correlation function subjected to the normalization constraint. From this minimization, we 



atedby 



16|,ll7|. 



get a set of Euler-Lagrange differential equations. The correlation functions are calcu 
solving these differential equations and then, the two-body energy, £'2, is computed 
This leads to the equations of state used in our present work. 



III. CALCULATIONS OF PROTONEUTRON STAR STRUCTURE PROPER- 
TIES 

In this section, we calculate the structure properties of neutron star in the different stages 
just after its formation by numerically integrating the TOV equation. The structure prop- 
erties of neutron star in these stages such as maximum mass, radius and temperature profile 
are calculated using the modern microscopic equations of state derived from constrained 
variational method employing the Argonne Vig potential pLd] as the inter-nucleonic interac- 
tion. At low densities, we also consider a hot dense matter model and use the equations of 
state obtained by Gondek et al. |3] and Strobel et al. 
stages. 



19|. Here is our results for different 



A. Lepton rich protoneutron star 



Since a neutron star, at the beginning of its lifetime (protoneutron star), is opaque with 
respect to neutrinos therefore, it contains a high lepton fraction (Y), 0.3 — 0.4. Furthermore, 
at this stage of neutron star, the entropy per baryon (s) is nearly constant throughout the 
star, 1 — 2kB. Since, we are also dealing with uncharged neutron star matter, the lepton 
and proton fractions should be equal {3, 3, 4, 5]. 

Our results for the pressure of protoneutron star matter as a function of mass density 
are presented in Fig. 1 at entropies s = Iks and s = 2kB with Y = 0.4 and Y = 0.3. It is 
seen that for a fixed value of entropy, the pressure of protoneutron star matter increases by 
decreasing the lepton fraction. We see that for a given value of lepton fraction, the equation 
of state for s = 2.0kB is stiffer than for s = LOks- In Fig . 1, We have compared our results 
with those of Gondek et al. [l^ and Strobel et al. 



19[. We have seen that the results of 



Gondek et al. [18| are nearly in agreement with our results for s = 2.0kB and Y = 0.3. But, 
the equations of state calculated by Strobel et al. [19] are stiffer than those of ours. 

In Fig. 2, we have presented the gravitational mass of protoneutron star as a function 
of central mass density for different values of entropy and lepton fraction. From Fig. 2, we 
can see that the gravitational mass increases by increasing the central mass density. It is 
shown that for different cases of entropy and lepton fraction, the gravitational mass exhibits 
different mass limits (maximum mass) as given in Table 1. From this table, the decreasing 
of maximum mass by decreasing the entropy and increasing the lepton fraction is seen. In 
Fig. 2 and Table 1, the results of Gondek et al. and Strobel et al. [loj] are also given 
for comparison. It is seen that these results are different from our calculations. 

The radius of protoneutron star as a function of central mass density is shown in Fig. 
3 for different entropies and lepton fractions. It is seen that as the central mass density 
increases, the radius decreases very rapidly and then reaches a nearly constant value. The 
radius of protoneutron star corresponding to the maximum mass for different entropies and 
lepton fractions is given in Table 2. It can be seen that the stiffer equation of state leads to 
the relatively higher radius. 

The variations of temperature throughout the protoneutron star matter with respect to 
the radial coordinate (r) are presented in Figs. 4 and 5 for s = Iks and s = 2kB with 
Y = 0.4 and Y = 0.3. We can see from the Figs. 4 and 5 that the temperature is decreasing 
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slowly by increasing r, but it drops rapidly near the protoneutron star crust. 



B. Beta-stable protoneutron star 

After complete deleptonization, neutrino trapped within the hot interior matter of neu- 
tron star do not affect the beta stability condition and therefore the lepton fraction is 
determined from the beta-equilibrium criteria 

mm 

. In this stage, we have calculated 
the structure properties of protoneutron star for both isentropic and isothermal paths. 

1. Isentropic paths 

For different values of entropy, our calculated equations of state of beta-stable protoneu- 
tron star are shown in Fig. 6. It is seen that the pressure of beta-stable protoneutron star 



matter increases by increasing entropy. In Fig. 6, the results of Strobel et al. 19|| are also 



plotted for comparison. There is a compatible difference between our results and those of 



Strobel et al. 



19|. 



Our calculated gravitational mass of beta-stable protoneutron star as a function of central 
mass density for s = LOks and s = 2.0kB is presented in Fig. 7. From this figure, the 
increasing of gravitational mass by increasing both central mass density and entropy is seen. 
As it is shown in Fig. 7, for higher values of central mass density, the gravitational mass 
shows a limiting value (maximum mass). The maximum mass of beta-stable protoneutron 
star is given in Table 3 for different entropies. It is seen that the maximum mass decreases 
by decreasing the ent.op^ In Fig. 7 and Table 3, we have compared on. .e=nlt= with the 
results of Strobel et al. p^]. It is seen that our results are different from those of Strobel et 



al. 

In Fig. 8, our results for the radius of beta-stable protoneutron star as a function of central 
mass density are plotted at different values of entropy. We see that the radius decreases very 
rapidly by increasing the central mass density and then exhibits a nearly constant value. The 
radius corresponding to the maximum mass of beta stable protoneutron star for s = LOfc^ 
and s = 2.0kB are also given in Table 3. This shows the decreasing of radius with respect 
to decreasing entropy. 

The temperature of beta-stable protoneutron star matter as a function of radial coordi- 
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nate (r) is shown in Fig. 9 for different entropies. The same as Fig.s 4 and 5, we see that 
the temperature decreases slowly by increasing r, but drops rapidly near the beta-stable 
protoneutron star crust. 

2. Isothermal paths 

We are now able to do the above calculations for the protoneutron star in the beta- 
equilibrium case at different isothermal paths. 

Our results for the equation of state, gravitational mass and radius of beta-stable pro- 
toneutron star at different values of temperature are presented in Figs. 10-12, respectively. 
At different temperatures, we have extracted the maximum mass and corresponding radius 
from these figures. Our results are given in Table 4. From this table, we can see that 
the maximum mass and radius of beta-stable protoneutron star decrease by decreasing the 
temperature. 

IV. SUMMARY AND CONCLUSION 

The protoneutron stars which are lepton rich and hot objects are formed just after the 
supernova explosion. In this paper, we have integrated TOV equation to compute the max- 
imum mass, radius and temperature profile for this object. As we have shown at a fixed 
value of lepton fraction, the mass and radius of protoneutron star decreases by decreasing 
entropy. The temperature is nearly constant in the core of protoneutron star and drops 
rapidly near its crust. The same properties of protoneutron star at beta equilibrium condi- 
tion have been also calculated along the isentropic and isothermal paths. We have shown 
that the maximum mass and corresponding radius increase by increasing temperature. 
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TABLE I: Our results for the maximum mass (Mq) of protoneutron star at different values of 
entropy (s) and lepton fraction (Y). The results of Gondek et al. [18] (GHZ) and Strobel et al. 
[19] (SSW) are also given for comparison. 





s = IMb 


s = 2.0kB 


Y = 0.4 


1.55 


1.56 


Y = 0.3 


1.63 


1.70 


SSW, Y = 0.4 


2.01 


2.03 


GHZ, Y = 0.4 




1.91 



TABLE II: Our results for the radius {km) of protoneutron star corresponding to the maximum 
mass at different entropies and lepton fractions. 





s = 1Mb 


s = 2Mb 


Y = 0.4 

Y = 0.3 


7.62 
8.13 


7.68 
8.35 



TABLE IV: Our results for the maximum mass {Mq ) of beta-stable protoneutron star and corre- 
sponding radius {km) at different temperatures. 





T 


= 5.0MeF 


T = 


lO.OMeV 


T = 


20.0MeF 


Mass 




2.19 




2.20 




2.23 


Radius 




7.85 




7.88 




8.01 



TABLE III: Our results for the maximum mass {Mq ) of beta-stable protoneutron star and cor- 
responding radius {km) at different entropies. The results of Strobel et al. [19] (SSW) for the 
maximum mass of beta-stable protoneutron star are also given for comparison. 





s = 1Mb 


s = 2Mb 


Mass 


1.68 


1.76 


Radius 


8.05 


8.15 


SSW 


1.98 


2.03 
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FIG. 1: Our results for the pressure {10^^ dyn/cm'^) of the protoneutron star matter versus mass 
density (lO^^^r/cm^) at entropies s = Iks and s = 2kB with lepton fractions Y = 0.4 and 0.3. The 
results of Gondek et al. [18] (GHZ) and Strobel et al. [19] (SSW) are also given for comparison. 




FIG. 2: As Fig. 1, but for the gravitational mass {Mq) of protoneutron star versus central mass 
density (lO^^gr/cm'^). 
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FIG. 3: Our results for the radius (km) of protoneutron star versus central mass density 
(W^^ gr / cm'^) at entropies s = Iks and s = 2kB with lepton fractions Y = 0.4 and 0.3. 
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FIG. 4: Our results for the temperature profile of protoneutron star at different entropies for 
Y = 0.4. 
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FIG. 5: As Fig. 4, but for Y = 0.3. 
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FIG. 6: Our results for the pressure {IQ^'^dyn/ cm?) of beta-stable protoneutron star matter as a 
function of mass density {10^'^ gr / cm?) at different entropies. The results of Strobel et al. [19] 
(SSW) are also given for comparison. 
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FIG. 7: As Fig. 6, but for the gravitational mass {Mq) of beta-stable protoneutron star versus 
central mass density (10^'^ gr / cm^) . 
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FIG. 8: Our results for the radius {km) of beta-stable protoneutron star versus central mass density 
(W^'^ gr / cm^) at different entropies. 
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FIG. 9: Our results for the temperature profile of beta-stable protoneutron star at different en- 
tropies. 
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FIG. 10: Our results for the pressure {10^'^ dyn / cm?) of beta-stable protoneutron star matter as a 
function of mass density {10^'^ gr / cm'^) at different temperatures. 
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FIG. 11: As Fig. 10, but for the gravitational mass {Mq) of beta-stable protoneutron star versus 
central mass density {W^'^gr /cm^). 
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FIG. 12: As Fig. 10, but for the radius (km) of beta-stable protoneutron star versus central mass 
density {10^'^ gr / cm'^) . 
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